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ABSTRACT 


A  collection  of  new  findings  are  presented  contributing  to  fundamental  physical  understanding  of 
broadband  (3-300  KHz)  transient  seismoacoustic  wave  phenomena  in  various  sand  conditions  (water- 
saturated,  naturally  deposited,  disturbed,  drained,  liquefied,  compacted,  and  water-covered)  important 
for  acoustic  modeling  of  seafloor  and  predicting  the  penetration  and  conversion  of  seismoacaoustic 
waves  in  sediments  The  experimental  results  compare  waveforms  of  different  types  of  seismoacoustic 
waves.  An  explanation  is  given  for  the  anomalous  slow  compressional  wave  observed  in  water-saturated 
sand  by  Boyle  and  Chotiros  (1992).  Other  experimental  results  presented  include  a)  the  effect  of 
ripples  on  surface  and  interface  waves  dispersion,  b)  the  existence  of  high-frequency  fast  compressional 
waves  (1650  m/s)  in  drained  sand,  c)  the  effect  of  puddles,  water  edges,  and  liquefied  regions  on 
converting  interface  waves  into  body  shear  waves,  f)  the  nonexistence  of  a  second  Rayleigh  wave  along 
with  Scholte  wave  in  soft  sediments,  and  g)  the  first  experimental  results  on  wedge  waves  along  the 
edge  of  sand  dollar  skeletons. 
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SUMMARY 


Little  is  known  about  the  seismoacoustic  properties  of  water-saturated  and  partially 
saturated  sandy  marine  sediments  in  the  frequency  range  1-100  KHz,  Sand  is  an  unstable 
medium  and  can  change  from  a  liquefied  state  to  a  varying  stiff  matrix  depending  on 
vibration,  water  content,  pressure,  air  bubbles,  and  aging.  In  the  surf  zone,  the  acoustic 
conditions  are  complicated  because  of  several ,  factors  (water  dynamics,  air  bubbles  in 
sediments  and  water,  sand  ripples,  sea  shells,  buried  shell  fragments,  tides,  density  and 
velocity  gradients,  sand  dollar  beds,  sorted  particles,  settling  of  transported  sediments, 
puddles  and  water  edges,  and  biological  activities)  as  depicted  in  Fig.  1. 

It  is  essential  to  study  the  fundamental  seismoacoustic  properties  of  water-saturated  and 
partially  saturated  sand  under  controlled  laboratory  settings  to  determine  the  different 
types  of  seismoacoustic  waves  existing  under  various  sand  conditions,  and  validate  current 
theoretical  models.  A  number  of  unresolved  controversial  problems  have  been  reported 
on  the  60  KHz  anomalous  slow  compressional  wave  (1200  m/s)  observed  by  Boyle  and 
Chotiros  [1]  in  water-saturated  sand,  and  the  acoustic  penetration  of  sandy  sediments  at 
incidence  angles  greater  than  Snell's  critical  angles  [1-2].  The  problem  is  complex  because 
a)  the  shear  wave  velocity  in  water-saturated  sand  increases  rapidly  with  depth  as  shown 
in  Fig.  2  indicating  that  penetrating  acoustic  waves  would  experience  strong  upward 
refraction' and  conversion,  b)  at  60  KHz,  the  shear  wavelength  is  about  5  grains  of  sand, 
and  the  validity  of  theoretical  transmission  coefficient  calculations  are  questionable  when 
the  grain  dimension  is  comparable  to  the  wavelength,  and  c)  the  seismoacoustic 
properties  of  sandy  sediments  are  not  understood  in  the  frequency  range  3-80  KHz. 
Another  current  controversial  fundamental  problem  is  the  physical  interpretation  of  the 


Fig.  1.  Diagram  depicting  examples  of  factors  affecting  seismoacoustic  properties  of 
sandy  seafloor  in  the  surf  zone. 
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Fig.  2.  Measured  shear  wave  velocity  at  different  depth  in  water-saturated  sand  using  a 
broadband  pulse  with  35  KHz  center  frequency. 

roots  of  the  dispersion  equation  corresponding  to  the  Rayleigh  and  Scholte  waves  [3-4]  at 
a  fluid-solid  interface  when  the  fluid  velocity  is  between  the  shear  wave  velocity  and  the 
compressional  wave  velocity  in  the  solid  (e.g.  sandy  sediments,  Plexiglas).  Only  a  Scholte 
wave  is  known  to  exist  at  such  interface,  however,  Padilla  et  al.  [3]  reported  theoretical 
and  experimental  results  apparently  indicating  the  coexistence  of  a  Rayleigh  wave  and  a 
Scholte  wave  at  a  water/Plexiglas  interface.  Several  new  seismoacoustic  findings  were 
obtained  while  investigating  the  controversial  problems  mentioned  above.  The  goal  of  the 
research  was  to  develop  better  fundamental  physical  understanding  of  seismoacoustic 
phenomena  in  sandy  sediments.  The  efforts  were  focused  on  identifying  gaps  in  current 
understanding  and  suggest  future  quantitative  research  needed  for  accurate  modeling  of 
seafloor  and  seismoacoustic  detection  of  buried  objects. 

The  new  findings  are  summarized  in  five  different  sections.  SECTION  1  provides  an 
explanation  [5]  for  the  time  of  arrival  of  the  anopialous  slow  compressional  wave  [1] 
based  on  the  depth-dependence  of  the  shear  wave  velocity.  Several  hypotheses  were 
reported  in  the  literature  as  summarized  in  [5-13],  Roughness  along  the  water/sand 
interface  can  convert  near-grazing  acoustic  waves  into  seismoacoustic  waves  in  the 
sediment  [11,  14-15],  The  scattered  field  from  roughness  propagates  laterally  in  the 
sediment  with  the  fast  compressional  wave  velocity  ~  1 700  m/s  and  consists  of  multiple 
delayed  components  spread  over  a  wide  time  window.  The  propagation  time  of  a  distinct 
1200  m/s  wave  could  not  be  explained  from  randomly  distributed  roughness.  Fig.  3 
demonstrates  the  effect  of  a  periodic  interface  roughness  on  generating  multiple  delayed 
fast  compressional  waves  using  a  laboratory  water/Plexiglas  ultrasonic  model.  The 
presence  of  surface  roughness  also  introduces  dispersion  of  interface  waves.  SECTION 
2  shows  the  effect  of  ripples  on  the  dispersion  of  Rayleigh  and  Scholte  waves.  The  results 
revealed  that  a  corrugated  surface  can  decrease  the  velocity  of  Rayleigh  and  Scholte 
waves  by  more  than  70%  [16  ]. 
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In  SECTION  3,  waveforms  from  various  sand  conditions  (water- saturated,  naturally 
deposited,  disturbed,  drained,  liquefied,  compacted,  and  water-covered)  are  compared. 
Two  surprising  results  were  obtained.  The  first  was  that  the  measured  velocity  of  high- 
frequency  (100  KHz)  compressional  waves  in  naturally  settled  water-saturated  sand 
remained  near  1700  m/s  when  the  sand  was  drained  and  did  not  drop  to  ~200  m/s  as 
anticipated  from  the  results  of  Velea  [17  ],  and  Tavossi  and  Tittman  [18  ]  near  10  KHz. 
The  decrease  in  density  as  air  bubbles  replaced  the  water  in  drained  sand  may  have  been 
counteracted  by  an  increase  in  stiffness  of  the  matrix  due  to  capillary  forces.  Further 
quantitative  research  in  needed.  The  second  surprising  result  was  the  coexistence  of  a 
broadband  fast  compressional  wave  and  a  slow  dispersive  low-frequency  compressional 
wave  in  water-covered  disturbed  sand  with  trapped  air  bubbles.  It  is  not  yet  understood 
whether  a  dispersive  Biot  wave  is  superimposed  over  this  slow  compressional  wave  in 
coarse  sand.  SECTION  4  shows  the  effect  of  water  edge  and  puddles  on  converting 
interface  waves  into  shear  waves.  Other  experimental  results  demonstrate  the  existence 
of  only  one  Scholte  wave  at  a  water/Plexiglas  interface  contradicting  the  results  of 
Padilla  et  al.  [3],  and  provided  an  explanation  for  the  excessive  amplitude  decay  of  the 
apparent  "Rayleigh"  wave  component  described  in  [3].  SECTION  5  presents  the  first 
experimental  results  on  nondispersive  flexural  wedge  waves  along  the  edge  of  sand  dollar 
skeletons  that  may  play  a  role  in  underwater  acoustic  interaction  with  a  seafloor  covered 
with  a  single-age  group  sand  dollar  bed.  In  the  CONCLUSION,  several  unanswered 
questions  are  presented  suggesting  needed  future  quantitative  research  leading  to  reliable 
detection  of  buried  objects. 


Fig.  3.  Effect  of  periodic  roughness  on  converting  near-grazing  underwater  acoustic 
waves  into  laterally  propagating  compressional  waves  in  elastic  bottom.  Top:  Received 
waveform  from  flat  interface.  Bottom:  Waveform  from  periodic  interface. 
(Water/Plexiglas  model  with  spatial  wavelength  6.34  cm  and  groove  depth  3  mm, 
compressional  velocity  =  2701  m/s,  range  =  15  cm,  horizontal  time  scale  10  ps/div.,  time 
at  beginning  of  trace  40  ps). 
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SECTION  1 


Anomalous  Slow  Compressional  Wave  And  Acoustic  Penetration 


1.1  An  Explanation  for  the  Anomalous  Slow  Wave  in  Water-Saturated  Sand 
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Letters _ 

An  Explanation  for  the  Anomalous  Ultrasonic  Slow  Wave  in 

Underwater  Sand 

Jacques  R.  Charauel,  Member,  IEEE 


Abstract — An  explanation  is  given  for  the  propagation 
time  of  the  well-known  anomalous  ultrasonic  slow  wave  ob¬ 
served  in  water-saturated  sand  using  a  three-layer  elastic 
model.  The  rapid  increase  of  elastic  properties  of  sand  with 
depth  causes  conversion  of  near-grazing  underwater  acous¬ 
tic  waves  into  multiple  coupled  shear  and  compressional 
waves. 


I.  Introduction 

Laboratory  experimental  studies  hv  Muir  et  al.  [1]  on 
the  transmission  of  near-grazing  underwater  ultrasonic 
waves  into  a  smooth  sandy  bottom  revealed  that,  unex¬ 
plained  acoustic  wave  peuel  nil  ion  occurred  at  angles  of  in¬ 
cidence  greater  than  the  critical  angle.  Chotiros  [2],  Boyle 
and  Chotiros  [3],  and  Altcnburg  and  Chotiros  [4]  observed 
an  anomalous  ultrasonic  wave  in  water-saturated  sand  that 
seemed  to  propagate  with  a  velocity  near  1200  m/s  at 
60  KHz.  Chotiros  [2]  suggested  that  the  anomalous  wave 
was  the  slow  compressional  wave  in  porous  media  pre¬ 
dicted  by  Biot  [5].  and  proposed  new  poroelastic  parameter 
values  because,  according  to  Biot's  theory,  the  slow  wave 
in  water-saturated  sand  should  be  highly  attenuated  and 
should  have  a  velocity  about  three  times  slower  than  this 
anomalous  wave.  Stoll  [7]  commented  on  inconsistencies 
related  to  the  input  parameters  to  the  Biot  theory  chosen 
by  Chotiros  [2],  The  Biot  slow  wave  was  not  conclusively 
detected  in  water-saturated  sand  [8]. 

Substantial  research  was  conducted  on  the  anomalous 
slow  acoustic  wave  and  on  the  poroelastic  properties  of  un¬ 
consolidated  sand  [2]— [4] .  [6]--[14]  (partial  list).  Researchers 
came  up  with  different  hypotheses  as  summarized  by  Ma- 
guer  et  al .  [10].  The  four  possible  explanations  given 
were:  poroelastic  Biot’s  theory  and  its  associated  second 
slow  wave,  water/sediment  interface  roughness  [11],  nar¬ 
row  beam  effects  [9],  and  evanescent  wave  scattering  by 
volume  inhomogeneities  [12].  In  addition.  Yamamoto  et.  al. 
[13]  reported  that  disturbances  in  the  seabed  for  instru¬ 
mentation  act  as  vertical  waveguides  that  can  be  mistaken 
as  a  false  Biot  slow  wave.  The  aim  of  this  letter  is  to  pro¬ 
pose  a  layered  elastic  model  as  a  hypothesis  to  explain  the 
travel  time  characteristics  of  the  apparent  slow  wave. 

Manuscript  received  April  IS,  1998:  accepted  July  1,  1998.  Work 
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II.  Hypothesis 

The  hypothesis  given  is  that  the  anomalous  ultrasonic 
slow  wave  [2]  is  the  result  of  multiple  conversion  of  near¬ 
grazing  underwater  acoustic  waves  incident  on  a  sandy 
bottom  with  depth-dependent  elastic  properties.  The  shear 
wave  velocity  of  underwater  unconsolidated  sand  increases 
rapidly  with  depth  [14]  in  the  upper  few  centimeters  form¬ 
ing  a  strong  velocity  gradient  layer  5  to  15  cm  thick.  An 
underwater  acoustic  wave  incident  at  angles  greater  than 
the  compressional  wave  critical  angle  is  partially  converted 
into  a  shear  wave  in  the  sandy  bottom.  The  presence  of  the 
large  shear  velocity  gradient  couples  the  shear  and  corn- 
pressional  waves  at  different  depths,  It.  is  known  that  the 
elastic  properties  of  sand  vary  continuously  with  depth. 
The  sand  will  be  represented  here  with  an  approximate 
layered  elastic  model  having  only  shear  wave  velocity  vari¬ 
ation  with  depth  and  a  constant  compressional  wave  veloc¬ 
ity.  Wave  refraction  depends  on  velocity  ratios,  and  wave 
conversion  efficiency  depends  on  both  density  ratios  and 
velocity  ratios.  The  purpose  of  this  preliminary  study  is  to 
focus  on  propagation  time  characteristics. 

III.  Hypothetical  Model  Parameters 

A  sound  source  is  located  in  water  at  distance  h,  above 
the  sand  interface.  A  receiver  is  buried  in  the  sand  at  depth 
d.  The  horizontal  distance  between  source  and  receiver  is  .r. 
The  sand  is  modeled  as  three  parallel  elastic  layers  over  a 
half-space.  The  compressional  wave  velocity  cp  is  1700  m/s 
in  all  the  sand  layers.  The  shear  wave  velocity  is  60,  140, 
190,  and  220  m/s  consecutively  starting  from  the  top  layer. 
The  respective  layer  thicknesses  are  4.5,  5.5,  6,  and  oc  cm. 
The  assumed  shear  wave  velocities  as  a  function  of  depth 
are  higher  than  that  reported  for  a  typical  sand  bed.  The 
shear  wave  velocity  was  not  given  directly  by  Boyle  and 
Chotiros  [3] .  They  mentioned  that  the  shear  wave  velocity 
of  sandy  sediments  similar  to  that  used  in  their  experi¬ 
mental  tank  is  about  150  m/s.  The  sediment  settled  for  18 
years  under  3  m  of  water  in  the  tank,  and  it  may  be  pos¬ 
sible  that  the  shear  wave  velocity  is  higher  than  150  m/s. 
The  preliminary  hypothetical  model  described  above  used 
shear  velocities  between  60  and  220  m/s  to  come  up  with 
one  simplified  layered  elastic  model  matching  the  propa¬ 
gation  time  characteristics  observed  at  different  receiver 
depths. 
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Horizontal  Distance  x  (m) 


Horizontal  Distance  x  (m) 

Fig.  1.  Propagation  time  plotted  as  function  of  x.  Open  squares: 
Calculated  results  from  a  three-layered  elastic  model,  (a)  d  =  0.23  m, 
and  (b)  d  —  0.35  m.  Solid  lines:  Experimental  data  by  Boyle  and 
Chotiros  [3]  replotted  from  their  Figs.  6  and  7  (d  =  0.33  m  and 
0.614  m).  Dashed  lines:  Simulated  numerical  results  by  Kibblewhite 
and  Wu  [6]  using  Chotiros’  modified  Biot  parameter  values  for  d  = 
0.33  m. 

IV.  Results  and  Discussions 

The  salient  features  of  the  experimental  propagation 
time  plots  by  Boyle  and  Chotiros  [3]  were  digitized  and 
replotted  to  facilitate  their  comparison  with  calculated  re¬ 
sults  from  the  proposed  hypothetical  model.  Four  distinct 
coexisting  travel  times  were  observed  in  the  experimen¬ 
tal  propagation  time  plots  (see  Fig.  1).  Corresponding  nu¬ 
merical  results  by  Kibblewhite  and  Wu  [6]  based  on  Biot 
parameters  modified  by  Chotiros  [2]  are  also  plotted  in 
Fig.  1(a).  The  Kibblewhite  and  Wu  [6]  explanation  re¬ 
sulted  in  two  rather  than  four  coexisting  travel  times.  The 
correct  hypothetical  model  should  explain  the  anomalous 
slow  acoustic  wave  [3]  and  predict  the  four  coexisting  travel 
times  in  the  experimental  data  [3]. 


The  values  of  source  height  h  and  receiver  depth  d  used 
in  the  calculations  were  determined  directly  from  the  ex¬ 
perimental  propagation  time  plots  in  [3]  and  not  from  the 
presumed  values  given  by  Boyle  and  Chotiros  [3]  because  of 
several  discrepancies  between  the  measured  and  expected 
propagation  time.  For  example,  the  value  of  h  was  deter¬ 
mined  from  the  arrival  time  t  of  the  wave  onset  and  not 
from  the  arrival  time  of  the  peak  (displayed  in  Fig.  4  in 
[3]).  Using  c  =  1489  m/s  and  t  =  302  /us  from  the  experi¬ 
mental  plot,  h  equals  1489  x  302  x  10-6  =  0.45  m  instead 
of  the  given  0.49  m  in  [3]. 

Calculated  results  from  the  hypothetical  layered  elastic 
model  are  presented  in  Figs.  1(a)  and  (b)  for  two  cases 
d  =  0.23  m  and  d  =  0.35  m  corresponding  to  the  experi¬ 
mental  results  (Figs.  6  and  7  in  [3]),  respectively,  where  the 
given  values  of  d  were  0.33  m  and  0.614  m.  The  two  receiver 
depths  used  in  the  calculations  were  determined  from  the 
arrival  time  t  —  458  /us  and  t  =  517  /us  of  the  wave  onset 
at  x  =  0  (Figs.  6  and  7  in  [3]),  h  —  0.45  m,  c  =  1489  m/s, 
and  cp  —  1700  m/s.  Note  that  the  receiver  depth  was  mis¬ 
printed  as  0.8  m  instead  of  0.08  m  (in  Fig.  5  of  [3]).  Propa¬ 
gation  time  for  selected  wave  components  were  calculated 
using  Snell's  law.  The  notation  used  is  as  follows:  subscript 
w  represents  waves  in  the  water,  and  no  subscript  means 
that  the  compressional  P  wave  or  the  shear  S  wave  com¬ 
ponents  remained  of  the  same  wave  type  (no  conversion) 
for  the  remainder  of  the  wavepath.  For  example,  PwS\SiP 
represents  a  P  wave  in  water  converting  into  shear  waves  in 
the  first  and  second  layers  and  propagating  the  remainder 
of  the  path  as  a  compressional  wave.  The  calculated  results 
were  superimposed  on  the  propagation  time  plots  by  Boyle 
and  Chotiros  [3].  The  four  sets  of  curves  in  Figs.  1(a)  and 
1(b),  from  bottom  to  top,  represent  the  converted  PWP. 
PwPiS^P-.  PWS\P,  and  PU,S\S2S$P  components,  respec¬ 
tively,  for  both  cases  d  —  0.23  m  and  d  —  0.35  m.  Cal¬ 
culated  propagation  time  results  obtained  from  a  single 
approximate  layered  model  simultaneously  matched  all  of 
the  experimental  multiple  propagation  time  features  pre¬ 
sented  in  [3]  for  different  source/receiver  geometries.  The 
preliminary  hypothetical  layered  model  presented  here  is 
not  unique;  it  was  presented  simply  to  demonstrate  that 
multiple  wave  conversion  also  can  provide  an  explanation 
for  the  anomalous  wave  without  contradicting  the  clas¬ 
sic  Biot  parameter  values  for  sand.  Other  more  elaborate 
five-layer  and  seven-layer  elastic  models  were  also  obtained 
by  the  author  matching  the  experimental  results  by  Boyle 
and  Chotiros  [3].  The  effect  of  continuous  velocity  gradient 
in  the  sand  and  the  relative  amplitudes  of  the  converted 
wave  components  in  underwater  sand  will  be  discussed  in 
a  separate  paper. 
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1.2  Experimental  Results  on  Refracted  Converted  Waves 

Ultrasonic  modeling  techniques  [19-22]  were  used  to  obtain  qualitative  and  quantitative 
experimental  results  under  controlled  laboratory  conditions.  An  experiment  was  first 
conducted  to  see  if  a  slow  Biot  wave  or  a  1200  m/s  "anomalous"  wave  could  be  detected 
in  water-saturated  sand.  The  experimental  model  consisted  of  a  10  cm  diameter  thin 
walled  plastic  tube  filled  with  water-saturated  sand.  Source  and  receiver  transducers  were 
placed  at  a  depth  of  43  cm  below  the  sand  surface.  The  source  was  excited  with  a 
broadband  pulse.  The  received  waveform  was  displayed  on  an  oscilloscope  with  a 
horizontal  time  scale  of  20  us/div.  The  results  showed  only  one  compressional  wave 
propagating  with  a  velocity  of  1652  m/s  as  shown  in  Fig.6,  The  waveform  did  not  reveal  a 
Biot  wave  or  the  "anomalous"  1200  m/s  slow  compressional  wave. 

The  effect  of  continuous  shear  wave  velocity  increase  with  depth  in  water-saturated  sand 
on  refracted  converted  waves  can  be  approximated  using  a  discretely  layered  model.  Fig. 
7  illustrates  refracted  converted  waves  at  a  single  solid/solid  interface  using  a 
Lexan/Plexiglas  experimental  model.  A  shear  wave  incident  from  the  Lexan  generated  a 
large  refracted  converted  compressional  wave  "SP"  in  the  Plexiglas.  The  model 
parameters  are  given  in  the  figure.  The  measured  elastic  properties  of  the  modeling 
materials  are  shown  in  parenthesis  in  the  following  order  density,  compressional  wave 
velocity,  and  shear  wave  velocity.  The  Lexan/Plexiglas  shear  wave  velocity  ratio  is 
879/1367.  Another  similar  experiment  was  conducted  on  a  wax/Lexan  interface  with  a 
smaller  shear  velocity  ratio  (600/879)  also  ^yielded  a  large  shear/compressional 
conversion  "SP"  as  shown  in  Fig.  8.  In  water-saturated  sand,  refracted  converted  waves 
are  expected  to  occur  at  a  much  lower  level.  Preliminary  refraction  conversion  results 
from  water-covered  sand  could  not  be  interpreted  from  the  known  fast  compressional 
wave  and  shear  wave  velocities.  The  puzzling  results  led  to  redirecting  the  research  to 
investigate  the  fundamental  seismoacoustic  properties  of  water-covered  sand  under 
various  conditions  as  described  in  SECTION  2. 
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Fig.  6.  Compressional  wave  in  water-saturated  sand. 
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SECTION  2 


2.  Seismoacoustic  Waves  in  Various  Sand  Conditions 

While  investigating  the  problem  of  refracted  converted  waves  in  water-covered  sand, 
unexplained  preliminary  results  were  obtained  that  influenced  the  course  of  the  research  to 
acquire  fundamental  understanding  of  seismoacoustic  properties  of  wet  sand.  These 
preliminary  refraction  conversion  results  could  not  be  interpreted  from  the  known  fast 
compressional  wave  and  shear  wave  velocities  (Fig.  9).  In  this  section,  measured 
waveforms  gathered  from  various  sand  conditions  are  presented. 

First,  here  is  a  brief  summary  of  the  findings  discussed  in  this  section.  Laboratory 
experimental  results  revealed  the  coexistence  of  a  fast  and  a  slow  compressional  wave  in 
water-covered  sand  under  certain  conditions.  The  slow  compressional  wave  had  a  time- 
varying  dispersion  and  its  velocity  initially  decreased  with  increasing  water  content,  and 
increased  with  repeated  stirring  and  compaction.  The  slow  compressional  wave  was 
more  detectable  in  disturbed  sand  and  in  coarse  sand.  It  was  surprising  to  observe  that 
the  velocity  of  the  fast  high-frequency  compressional  wave  in  water-saturated  sand 
remained  practically  unchanged  when  the  sand  was  drained.  This  indicated  that  the  high- 
frequency  components  propagated  through  the  sand  grain  matrix  and  were  little  affected 
by  the  drained  water  that  existed  in  the  pores.  Current  theories  may  need  to  be  modified 
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Fig.  9.  Preliminary  results  showing  waveforms  in  wet  sand  recorded  at  a  fixed  horizontal 
range  (5.5  cm)  and  different  receiver  depth  (0.25  cm,  3.25  cm,  and  5,25  cm).  Shear 
source  and  compressional  receiver.  Bottom  trace  shows  dispersion  and  an  unexplained 
refracted  wave.  Horizontal  time  scale  200  ys/div. 
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to  account  for  this  effect.  Complete  wetting  of  sand  grains  is  not  always  achieved 
immediately  upon  immersing  the  sand  in  water  resulting  in  unstable  time-varying 
seismoacoustic  properties.  The  new  findings  raised  several  questions  that  are  presented  in 
the  CONCLUSION  section. 

Examples  are  given  in  Fig.  10  demonstrating  typical  variations  in  the  received  waveform 
from  water-covered  sand  under  various  conditions.  Fig.  10(a)  shows  a  fast  high- 
frequency  compressional  wave  (1650  m/s)  followed  by  a  low-frequency  shear  wave  (140 
m/s)  detected  from  naturally  settled  compacted  sand.  The  waveforms  in  Fig.  10(b)-(e) 
were  obtained  from  water-covered  sand  that  was  previously  disturbed  while  the  water 
level  was  below  the  sand.  The  high-frequency  fast  compressional  wave  was  highly 
attenuated,  and  the  shear  wave  remained  practically  unchanged.  A  new  slow 
compressional  wave  (190-400  m/s)  appeared  in  the  disturbed  sand  arriving  before  the 
shear  wave.  Its  amplitude  and  arrival  time  varied  greatly  depending  on  the  effect  of 
disturbance  and  trapped  microbubbles.  The  source-receiver  distance  was  10  cm,  and  the 
source  was  excited  with  a  single  broadband  pulse. 


Fig.  10.  Waveform  variations  in  water-covered  sand:  (a)  naturally  deposited  sand,  (b)- 
(e)  water-covered  sand  disturbed  when  water  level  was  below  sand.  Notice  presence  of 
a  high-frequency  fast  compressional  wave  1650  m/s  in  (a)  followed  by  a  low-frequency 
shear  wave  140  m/s.  The  fast  compressional  wave  was  highly  attenuated  in  disturbed 
sand  while  the  shear  wave  was  practically  not  affected.  A  new  slow  compressional  wave 
appeared  in  water-covered  disturbed  sand  190-400  m/s.  Source-receiver  distance  was  4 
cm.  Time  scale  50  pis/div.. 
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Fig.  11  compares  the  waveforms  of  compressional  waves  in  drained  sand  and  water- 
covered  sand.  The  compressional  wave  velocity  in  drained  sand  was  190  m/s  with  a  10 
KHz  center  frequency.  The  waveform  from  the  water-covered  sand  exhibited  two 
compressional  wave  components  consisting  of  a  high-frequency  fast  compressional  wave 
(~  1650  m/s)  and  a  slow  low-frequency  compressional  wave.  The  shear  wave  velocity  in 
drained  sand  was  greater  than  the  shear  wave  in  water-covered  sand  (Fig.  12).  Variations 
in  the  shear  wave  amplitude  and  arrival  time  due  to  draining,  saturation,  liquefaction,  and 
compaction  are  compared  in  Fig.  12(a)-(d)..  The  shear  wave  velocity  decreased  with 
liquefaction  and  water  loading.  The  horizontal  time  scale  was  200  ps/div. 

Waveforms  of  the  fast  high-frequency  compressional  wave  are  examined  in  Fig.  13  using  a 
5  ps/div.  time  scale.  The  fast  compressional  wave  velocity  in  water-covered  naturally 
deposited  sand  was  1610  m/s  (Fig,  13(a)).  The  velocity  dropped  from  1610  m/s  to  1556 
m/s  with  liquefaction  Fig.  13(c).  One  surprising  result  was  obtained  when  the  water  level 
was  dropped  9  cm  below  the  sand  and  the  fast  compressional  wave  did  not  change 
significantly  (Fig.  13(b)).  A  similar  experiment  was  conducted  using  sea  water  and  it 
yielded  the  same  surprising  result  (Fig.  14).  The  measured  densities  of  the  dry  ,  water- 
saturated,  and  drained  sand  were  1765,  2169,  and  2148  kg/m3 ,  respectively. 


Drained  Wet  Sand 


Water-covered 


Fig.  1 1  .  Comparison  of  compressional  waves  in  drained  wet  sand  and  water-covered 
sand.  The  top  trace  shows  a  slow  -10  Khz  compressional  wave  (190  m/s),  The  bottom 
trace  shows  two  compressional  waves  one  fast  high  frequency  compressional  wave  and 
one  slow  low-frequency  wave  (faster  than  the  shear  wave).  The  source  was  excited  with 
one  broadband  pulse.  Time  scale  100  ps/div.. 
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Fig.  12.  Variations  of  shear  wave  amplitude  and  arrival  time  in  naturally  deposited 
underwater  sand  due  to  liquefaction  and  compactness  The  high-frequency  components  at 
the  beginning  of  the  traces  are  due  to  compressional  waves  Notice  shear  wave  velocity 
decreased  with  liquefaction  and  water  loading  Compaction  increased  the  shear  wave 
velocity  of  the  previously  liquefied  sand.  Source-receiver  distance  4.7  cm  at  2  cm  depth. 
Horizontal  time  scale  200  ps/div., 


(a)  Water-Covered 
Naturally  Deposited  Sand 

(1610  m/s) 

(b)  Dropped  Water  Level  9  cm  Below  Sand 


(c)  Water-Covered  Liquefied  (1556  m/s) 


Fig.  13.  Comparison  of  high-frequency  compressional  wave  in  sand.  Compressional 
wave  velocity  dropped  from  1610  m/s  to  1  556  m/s  with  liquefaction.  Compressional  wave 
attenuation  decreased  in  liquefied  sand  Horizontal  time  scale  5  ps/div. 
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(a)  Water-saturated 


(b)  Drained 


Fig.  14.  Surprising  results  showing  the  fast  high-frequency  compressional  wave  velocity 
(1750  m/s)  in  water-saturated  sand  (sea  water)  remained  practically  unchanged  (dropped 
to  only  1676  m/s)  when  the  sea  water  was  drained  and  the  sand  was  sponge-dried. 
Source-receiver  horizontal  separation  distance  was  4  cm. 


The  center  frequency  of  the  broadband  pulse  used  to  excite  the  source  was  selected  to  enhance 
low  or  high  frequency  components  in  different  experiments.  The  fast  compressional  wave  in 
water-saturated  sand  (Fig.  15(b))  is  compared  with  the  compressional  wave  in  water  (Fig. 
15(a))  using  a  400  KHz  pulse  center  frequency.  The  propagation  distance  was  3.3  cm. 

The  high-frequency  fast  compressional  wave  results  presented  in  Fig.  13(b)  and  Fig.  14(b) 
were  surprising  because  the  compressional  wave  velocity  in  drained  sand  should  be  about  8x 
slower  than  the  velocity  in  saturated  sand  [17-18]  measured  near  10  KHz.  In  order  to  reverify 
the  results,  another  controlled  experiment  was  conducted  to  measure  the  fast  compressional 
wave  in  a  small  wet  sand  sample  1  x  2  x  4  cm  as  shown  in  Fig.  16.  The  procedure  was  as 
follows:  source  and  receiver  transducers  were  positioned  in  the  middle  of  a  water  tank.  The 
distance  between  the  transducers  was  fixed  at  4  cm.  Sand  was  sprinkled  very  slowly  over  the 
water  to  form  naturally  settled  sand.  The  tank  wall  was  tapped  periodically  to  compact  the 
sand  untill  the  tank  was  filled  with  water-saturated  sand.  The  source  was  excited  wih  one 
broadband  pulse.  The  received  fast  compressional  wave  was  captured  on  the  oscilloscope. 
The  excess  water  on  the  sand  surface  was  removed,  and  the  sand  was  carved  out  from  all  the 
tank  except  for  a  rectangular  sample  1  x  2  x  4  cm  naturally  held  between  the  parallel  faces  of 
the  transducers.  During  the  carving  process,  the.  received  fast  compressional  wave  was 
monitored  and  it  did  not  change.  The  sample  was  resaturated  and  the  compressional  waveform 
was  recorded  as  shown  in  the  top  trace  of  Fig,  16.  The  measured  compressional  wave  velocity 
was  1675  m/s.  The  second  trace  was  obtained  after  sponge  drying  the  sand  sample.  Notice 
that  the  compressional  wave  time  of  arrival  remained  almost  unchanged.  The  sand  sample  was 
allowed  to  dry  using  blowing  air  on  the  sample.  As  the  sample  was  drying,  the  received 
waveform  was  recorded  periodically.  The  bottom  trace  in  Fig,  16  was  obtained  from  the 
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"almost  dry"  sample  as  the  sand  bar  crumbled.  This  experiment  confirmed  that  the 
compressional  wave  velocity  of  high-frequency  components  above  80  KHz  remained  close  to 
the  water-saturated  velocity  and  did  not  decrease  by  a  factor  of  eight.  It  seems  that  the  high- 
frequency  waves  propagated  in  the  naturally  settled  compacted  sand  as  in  a  consolidated  matrix 
not  significantly  affected  by  the  drained  water.  Theoretical  models  and  selected  Biot 
parameters  need  to  account  for  the  compressional  wave  behavior  decribed  above  [7-8, 25], 

Bachrach  and  Nur  [25]  calculated  the  compressional  wave  velocity  as  function  of  water 
saturation  in  sand  based  on  the  Biot-Gassman  prediction.  The  calculated  results  show  that  the 
compressional  wave  velocity  initially  decreases  as  the  water  content  is  increased,  and  then  it 
starts  to  increase  above  85%  water  saturation.  At  99%  water  saturation,  the  compressional 
wave  velocity  is  still  less  than  200  m/s,  According  to  the  theoretical  model,  the  last  1%  of 
water  saturation  increases  the  compressional  wave  velocity  by  almost  a  factor  of  8. 


(a)  Water 


(b)  Saturated  Sand 


Fig.  15.  Received  compressional  wave  in  water  (top  trace),  and  in  water-saturated  sand 
(bottom  trace).  The  distance  between  source  and  receiver  was  3.3  cm. 
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SOURCE 
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■maw 


ALMOST  DRY 


Fig.  16.  Effect  of  decreasing  water  content  in  wet  sand  on  compressional  wave  .  Waves 
propagated  and  detected  in  a  wet  sand  block  1x2x4  cm  carefully  carved  from  naturally 
deposited  sand.  The  compressional  wave  velocity  remained  almost  constant  as  the  water 
content  decreased  over  a  wide  range  of  moisture  content.  Source-receiver  distance  4  cm. 
Horizontal  time  scale  10  ps/div. 
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The  arrival  time  of  the  slow  compressional  wave  relative  to  the  shear  wave  is  shown  in  Fig.  17, 
where  shear  source  and  receiver  transducers  were  utilized.  Fig.  17(a)  shows  the  slow 
congressional  wave  arriving  before  the  shear  wave  in  sprinkled  compacted  sand.  When  the 
sand  was  liquefied,  the  slow  compressional  wave  amplitude  decreased  and  new  high-frequency 
components  appeared  in  Fig.  17(b).  Experiments  using  different  types  of  sand  were 
conducted.  The  slow  compressional  wave  was  observed  more  easily  in  coarse  sand  than  in  fine 
sand.  Waveforms  from  coarse  and  fine  sands  are  compared  in  Fig.  18. 
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Examples  on  the  effect  of  repeated  liquefaction  and  compaction  on  compressional  waves  in 
water-covered  sand  are  given  in  Fig.  19  and  Fig.  20.  The  source  pulse  waveform  is  displayed 
on  the  left  portion  of  Fig.  19.  With  successive  liquefaction  and  compaction,  the  amplitude  of 
the  high-frequency  fast  compressional  wave  increased,  and  the  velocity  of  the  dispersive  slow 
compressional  wave  also  increased  (trace  sequence  from  top  to  bottom).  The  bottom 
waveform  in  Fig.  20  shows  how  the  slow  low-frequency  compressional  wave  overlapped  a 
portion  of  the  fast  compressional  wave  as  its  velocity  increased  with  repeated  liquefaction  and 
compaction. 

It  was  also  observed  that  the  characteristics  of  compressional  waves  in  water-covered  sand 
vary  with  time  (Fig.  21  and  Fig.  22).  Fig.  21  compares  the  waveforms  at  5  minutes  and  65 
minutes  following  settling  and  compaction  of  the  sand.  The  velocity  of  the  slow  compressional 
wave  decreased  with  time.  Stirring  and  compacting  the  sand  increased  the  velocity  of  the  slow 
compressional  wave  as  demonstrated  in  Fig.  21(c). 

The  appearance  of  the  slow  compressional  wave  in  water-covered  sand  sometimes  occurs 
after  waiting  several  seconds  (~  30  seconds)  as  shown  in  Fig.23.  The  exact  cause  has  not  yet 
been  determined,  however,  it  may  be  due  to  the  formation  of  bubbles,  settling  of  fine  particles, 
or  changes  in  grain  to  grain  contact. 

The  observed  slow  compressional  wave  in  water-covered  sand  is  dispersive.  Measured 
dispersion  results  are  presented  in  Fig.  24  for  play  sand,  and  in  Fig.  25  for  silica  sand  (150  pim- 
425  pirn).  The  measured  dispersion  characteristics  are  time  dependent.  The  results  shown  in 
Fig.  25  were  obtained  30  minutes  after  settling  and  compaction  of  the  sand.  Quantitative 
studies  are  needed  to  characterize  the  spatial  distribution  of  water  content,  grain  to  grain 
wetting  contact  conditions,  and  microbubbles  content.  The  Biot  wave  has  not  yet  been 
identified  in  water-saturated  sand.  Hickey  and  Sabatier  [31]  reported  on  the  dispersion  of  a 
slow  compressional  wave  in  air-filled  sand. 


Fig.  19.  Effect  of  repeated  liquefaction  and  compaction  on  compressional  waves  in  water- 
covered  sand  (  x  =  3  cm),  The  source  broadband  pulse  waveform  is  shown  on  the  left. 
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Fig.  20.  Another  example  showing  effect  of  repeated  liquefaction  and  compaction  on 
congressional  waves  in  water-covered  sand  (x  =  3  cm).  Time  scale  20  ps/div.. 
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(a)  t  =  5  min. 


(b)  t  =  65  min. 

(c)  t  ~  65  min. 

+  Stirred  and  Compacted 


Fig.  21.  Effect  of  wetting  time  and  disturbance  on  compressional  waves  in  water-covered 
silica  sand  (425  pm-850  pm).  The  dry  sand  was  "rained"  (sprinkled)  over  the  water  tank 
and  was  compacted.  Time  measurement  began  immediately  after  settling.  Source-receiver 
range  was  4  cm.  Time  scale  50  ps/div.. 


Fig.  22.  Aging  effect  on  compressional  waves  in  water-covered  silica  sand  (425  pm-850 
pm).  Source  -receiver  distance  was  6  cm.  Horizontal  time  scale  100  ps/div.. 
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Fig.  23.  Appearance  of  the  slow  compressional  wave  occurs  sometimes  after  a  waiting 
time  of  several  seconds  (~  30  seconds)  in  water-covered  sand,  may  be  due  to  the 
formation  of  bubbles  or  settling  of  the  fine  particles.  Source-receiver  distance  was  2.96 
cm.  Horizontal  time  scale  50  ys/div.. 


Sonoquest  JRC-74-99 


25 


J.  R.  Ctuunuel 


FREQUENCY  (KHz) 

Fig.  24.  Plot  of  measured  dispersion  of  compressional  wave  phase  velocity  in  water- 
covered  play  sand  after  substantial  stirring  and  compaction.  Vertical  wave  propagation 
path  in  sand. 


FREQUENCY  (KHz) 


Fig.  25.  Plot  of  measured  dispersion  of  compressional  wave  phase  velocity  in  water- 
covered  silica  sand  (150  pm  -425  pm)  30  minutes  after  settling  and  compaction.  Vertical 
wave  propagation  path  in  sand. 
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SECTION  3 


3.1  Effect  of  Water  Level  on  Rayleigh  and  Shear  Waves 

Near  the  beach,  variations  in  water  level  greatiy  affect  the  Rayleigh  and  shear  waves. 
Laboratory  experimental  results  on  sand  are  presented  in  Fig.  demonstrating  variations  in 
received  waveforms  as  the  water  level  changed.  A  piezoelectric  source  transducer  was  placed 
on  the  flat  sand  surface  of  drained  wet  sand  to  generate  Rayleigh  waves.  A  2.5  cm  diameter 
shear  receiver  transducer  was  buried  at  a  distance  of  10  cm  flush  with  the  sand  surface  and 
oriented  horizontally  to  detect  Rayleigh  and  shear  waves.  The  Rayleigh  wave  in  drained  wet 
sand  is  shown  in  the  top  trace  having  a  velocity  of  65  m/s.  The  shear  wave  arrived  a  little 
before  the  Rayleigh  wave  and  had  a  higher  center  frequency. 

When  the  water  level  was  increased  to  reach  the  sand  surface,  both  the  shear  and  Rayleigh 
wave  amplitude  and  velocities  decreased.  The  measured  Rayleigh  wave  velocity  in  water- 
saturated  sand  was  44  m/s.  Two  other  waveforms  are  presented  in  Fig.  26  with  the  water  level 
at  1 .5  cm  and  3  cm  below  the  sand  surface.  The  ratio  of  shear  wave  velocity  to  Rayleigh  wave 
velocity  varied  nonlinearly  as  the  water  level  dropped.  Notice  as  the  water  level  dropped  3  cm 
below  the  sand  surface,  the  Rayleigh  wave  amplitude  (right  peak)  increased  and  its  propagation 
time  decreased.  The  wave  velocities  should  increase  with  depth  in  sand,  however,  the  presence 
of  the  water  level  at  a  short  distance  below  the  surface  decreases  the  velocity  of  the  lower 
frequency  components  reaching  the  deeper  water-saturated  sand.  Future  studies  are  needed  to 
quantify  these  water  level  effects. 


Fig.  26.  Effect  of  water  level  in  sand  on  Rayleigh  and  shear  waves. 
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3.2  Conversion  of  Interface  Waves  Into  Shear  Waves  by  Water  Edge 


The  conversion  of  interface  waves  into  shear  waves  by  a  water  edge  is  demonstrated.  A 
water/Plexiglas  laboratory  ultrasonic  model  was  used  in  the  experiment.  A  2  mm  diameter 
piezoelectric  transducer  was  placed  on  the  flat  suface  of  a  Plexiglas  block  simulating  an  elastic 
half-space.  A  shear  receiver  was  positioned  2  cm  below  the  surface  at  a  distance  of  25.7  cm. 
Fig.  27  shows  two  waveforms  obtained  from  the  Plexiglas  block  partially  immersed  in  water. 
The  angle  of  intersection  between  the  water  and  Plexiglas  varied.  The  top  waveform  shows 
the  direct  shear  wave  "SI"  followed  by  an  unexpected  large  shear  wave  "S2"generated  from 
partial  conversion  of  the  incident  Rayleigh  wave  by  the  water  edge.  The  angle  of  intersection 
was  90  degrees  (ignoring  the  curved  surface  tension  region).  The  distance  between  the  source 
and  the  water  edge  was  18.7  cm.  The  bottom  waveform  was  recorded  with  the  block  tilted  to 
make  a  10  degree  angle  of  intersection,  Again,  a  significant  shear  wave  was  generated  by  the 
conversion  of  the  Rayleigh  wave.  Similar  results  were  observed  from  the  conversion  of  a 
Scholte  wave  into  a  shear  wave  by  the  water  edge  when  the  source  was  located  underwater  on 
the  Plexiglas  surface.  The  results  hint  the  potential  role  of  liquefied  localized  regions  in 
seafloor  on  converting  interface  waves  into  shear  waves  contributing  to  transmission  loss  in 
very  shallow  water  even  in  the  absence  of  topography. 


Source 


Angle 


90  degrees 


1 0  degrees 


Fig.  27.  Shear  wave  S2  generated  in  elastic  bottom  from  the  conversion  of  a  Rayleigh  wave 
interacting  with  the  water  edge.  Horizontal  time  scale  10  us/div.  Time  at  beginning  of  trace 
160  us.  SI  is  the  direct  shear  wave  from  source  to  receiver. 
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The  presence  of  a  water  puddles  on  the  surface  of  a  solid  can  also  partially  convert  an  incident 
Rayleigh  wave  into  shear  waves.  Fig.  28  shows  the  effect  of  varying  the  diameter  of  a  single 
circular  water  puddle  placed  on  a  Plexiglas  surface  on  converting  an  incident  Rayleigh  into  a 
shear  wave.  The  waveform  was  recorded  with  the  puddle  diameter  0,  5,  10,  and  20  mm. 
Distance  from  source  to  puddle  18.7  cm.  horizontal  distance  from  puddle  to  receiver  7  cm, 
receiver  depth  2  cm. 


Source 


Water  Puddle 


/ 


Air 


Plexiglas 


Receiver 


Converted 
Shear  wave 


Fig.  28.  Effect  diameter  of  water  puddle  placed  on  solid  flat  surface  (Plexiglas)  on  converting 
incident  Rayleigh  waves  into  shear  waves  penerating  the  solid. 
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3.3  On  Coexistence  of  Rayleigh  and  Scholte  Waves 


A  current  controversial  fundamental  problem  needs  to  be  resolved  regarding  the 
coexistence  of  a  Rayleigh  wave  and  a  Scholte  wave  at  a  fluid/solid  interface  when  the  fluid 
sound  velocity  is  between  the  shear  and  compressiona!  wave  velocities  in  the  solid  [3] 
(e.g.  sandy  sediments,  Plexiglas,  sandstone).  Padilla  et  al.  [3]  presented  numerical  results 
showing  there  might  exist  two  surface  waves  from  the  roots  of  the  dispersion  equation  and 
confirmed  the  existence  of  a  Rayleigh  wave  and  a  Scholte  wave  from  water/Plexiglas 
experimental  results.  In  1982,  Chamuel  [38]  presented  experimental  results  from  an 
ultrasonic  water/Plexiglas  model  indicating  that  only  one  interface  Scholte  wave  exists  at 
such  interface  using  point  source  and  receiver  transducers  placed  on  the  interface. 

In  view  of  the  findings  presented  in  the  previous  section  on  the  generation  of  shear  waves 
from  the  conversion  of  Rayleigh  waves  by  a  water  edge,  the  interpretation  of  the 
experimental  results  by  Padilla  et  al.  [3]  from  water/Plexiglas  became  questionable. 

A  new  experiment  was  conducted  using  a  water/Plexiglas  model  that  provided  an 
explanation  for  the  excessive  amplitude  decay  of  the  apparent  "Rayleigh"  wave  component 
observed  by  Padilla  et  al.  [3],  The  new  findings  indicated  that  the  experimental  results  in 
[3]  were  probably  misinterpreted  by  referring  to  the  converted  shear  wave  as  a  Rayleigh 
wave  coexisting  with  the  Scholte  wave.  The  Rayleigh  wave  from  the  free  Plexiglas  model 
in  air  is  shown  in  the  top  trace  of  Fig.  29  (no  separate  shear  wave  component  observed). 
A  piezoelectric  transducer  (2  mm  diameter)  was  placed  on  the  surface  of  the  Plexiglas  to 
generate  a  Rayleigh  wave.  A  shear  receiver  transducer  was  coupled  to  the  end  of  the 
Plexiglas  block  .  The  distance  between  source  and  receiver  was  14  cm.  The  middle  trace 
shows  two  wave  components  (a  shear  wave  and  a  Scholte  wave)  generated  from  the 
conversion  of  the  Rayleigh  wave  by  partial  immersion  in  water  wetting  5  cm  length  of  the 
wavepath.  Touching  the  surface  of  the  Plexiglas  attenuated  the  Scholte  wave  and  did  not 
disturb  the  shear  wave  generated  from  the  water  edge  conversion  as  shown  in  the  bottom 
trace  verifying  that  the  remaining  unattenuated  component  was  a  body  shear  wave  and  not 
a  surface  wave.  Spherical  spreading  of  the  shear  wave  accounts  for  the  unexplained 
excessive  amplitude  decay  of  the  apparent  "Rayleigh"  wave  component  observed  by 
Padilla  et  al.  [3],  The  arrival  time  of  the  converted  shear  wave  by  the  water  edge 
decreased  as  the  immersion  length  increased. 

Fig.  30  shows  the  calculated  arrival  time  of  the  Scholte  wave  and  shear  wave 
components  generated  by  the  conversion  of  the  Rayleigh  wave  by  the  water  edge.  In  the 
figure,  the  experimental  results  by  Padilla  et  al. [3]  are  also  replotted  with  the  calculated 
arrival  time  of  the  hypothetical  second  Rayleigh  wave.  Notice  the  good  match  between 
the  experimental  results  from  Reference  [3]  and  the  new  calculated  results  based  on  the 
converted  shear  wave.  The  calculated  results  were  obtained  using  the  following 
parameters:  compressional  wave  velocity  2745  m/s,  shear  wave  velocity  1363  m/s, 
Rayleigh  wave  velocity  1295  m/s,  Scholte  wave  velocity  1058  m/s,  range  14  cm,  and 
receiver  depth  3  mm  (effective  depth  of  shear  receiver  being  on  the  side  of  the  Plexiglas 
block).  The  findings  are  important  for  interpreting  the  roots  of  the  dispersion  equation  and 
predicting  energy  partitioning  in  seismoacoustic  signals  in  sandy  sediments. 
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component  is  a  converted  shear  wave 


Fig.  29.  Experimental  results  demonstrating  the  conversion  of  a  Rayleigh  wave  by  partial 
immersion  of  a  flat  solid  in  water  produces  a  Schote  wave  and  a  body  shear  wave  and  not 
a  second  interface  wave  as  reported  in  [3], 


Fig.  30.  Calculated  arrival  time  of  Scholte  and  shear  waves  generated  from  the  conversion 
of  a  Rayleigh  wave  as  a  function  of  immersion  length  for  a  water/Plexiglas  model.  The 
experimental  results  by  Padilla  et  al.  [3]  are  plotted  (open  squares)  together  with  the 
calculated  arrival  times  of  the  hypothetical  second  Rayleigh  wave  (thin  solid  line). 
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SECTION  4.  Ripples 


4.1  Scholte  Wave  Dispersion  by  Rippled  Liquid/Solid  Interface  Topography 
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Scholte  Wave  Dispersion  by  Rippled  Liquid/Solid  Interface  Topography 

Jacques  R.  Chamuel 

Sonoquest  Advanced  Ultrasonics  Research,  P.  0.  Box  8U53,  Wellesley  Hills,  MA  02181-0001 

ict:  Quantitative  experimental  result*  are  presented  characterizing  the  dispersion  of  Scholte  wave*  propagating  aioog 
led  surface  of  an  immersed  "soft*  solid  half-space.  One  model  has  a  sinusoidal  profile  with  a  ratio  of  corrugation 
ude  to  wavelength  of  0,25.  Scholte  and  Rayleigh  wave  experimental  dispersion  results  are  compared  with  numerical 
i  on  Rayleigh  wave  dispersion  along  a  "hard'  corrugated  solid  half-space  (N.E.  Glass,  R.  Loudon,  and  A.  A. 
Judin,  Physical  Review  0  24,  6843-6861  (1981)),  The  findingi  reveal  that  a  rippled  "soft"  liquid/solid  interface  can 
asethc  velocity  of  high-frequency  Scholte  wave  components  propagating  normal  to  the  ripples  by  more  than  70ft. 


INTRODUCTION 

t  effect  of  a  rippled  liquid/solid  interface  on  the  dispersion  of  Scholte  waves  is  not  known.  Frequency- 
ndenl  multiple  scattering  phenomena  occur  as  near-grazing  underwater  acoustic  waves  and  Scholte  waves 
act  with  periodic  seafloor  topography  having  dimensions  comparable  to  the  wavelength.  In  order  to  develop 
,icai  understanding  of  the  interaction  of  broadband  transient  underwater  acoustic  waves  and  Scholie  waves  with 
>pled  ocean  bottom,  experimental  ultrasonic  experiments  were  conducted  on  immersed  isotropic  'soft'  models 
t  a  sinusoidally  corrugated  surface.  Quantitative  experimental  results  were  obtained  characterizing  the  dispersion 
Rayleigh  and  Scholte  waves  caused  by  a  periodic  tiquid/solid  interface. 


MODEL  EXPERIMENTS  AND  RESULTS 


A  ainuoidal  partem  was  machined  along  the  edge  of  a  4.5  mm  thick  Plexiglas  plale.  The  spatial  wavelength  of 
e  periodic  pattern  was  2  cm  and  the  peak-to-peak  amplitude  was  1  cm.  The  measured  Plexiglas  velocities  were 
;  follows:  compressions!  plale  velocity  =  23  J  3  m/s,  shear  wave  velocity  *  1344  m/s,  flat  edge  Rayleigh  wave 
elocity  *>  1220  m/s,  flat  edge  Scholte  wave  velocity  =  1109  m/s.  The  top  trace  of  Figure  1(a)  compares  the 
ispersed  waveforms  of  Rayleigh  and  Scholte  waves  propagating  along  the  periodic  surface.  Rayleigh  and  Scholte 
vaves  were  obtained  from  air/Plexiglas  and  water/Plexiglas  models,  respectively.  The  Scholte  waveform  was 
,-tretched  more  to  the  right  than  the  Rayleigh  wave  indicating  that  it  was  slower  than  the  Rayleigh  wave.  The 
oottom  trace  of  Figure  1(a)  shows  the  arrival  time  of  the  corresponding  Rayleigh  wave  along  a  flat  surface  (no 


(a)  (b)  (c) 


FIGURE  1.  Comparison  of  received  waveforms  and  arrival  lime;  (a)  Top  Truce;  Rayleigh  and  Scholte  waves 
along  corrugated  surface  (range  x  =  12  cm),  BoUom  Truce:  reference  Rayleigh  wave  from  flat  solid  surface,  (b) 
Rayleigh  wave  along  corrugated  surface  detected  at  two  positions  separated  by  one  corrugation  cycle  (2  cm),  (c) 
Scholte  wave  along  corrugated  surface  detected  at  the  same  two  positions  in  (b). 
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FREQUENCY  (XHi) 

trace  1  Theoretical  (Glass  et  al.  19B1)  Rayleigh  Wave 

'  *"  trace  2  Experimental  [Chamuel  1997)  Rayleigh  Wave 

trace  3  Experimental  [Chamuel  1997]  Scholte  Wave 

FIGURE  2.  Rayleigh  and  Scholte  wave  dispersion  experimental  results  compared  with  approximately  scaled 
theoretical  dispersion  results  by  Class  et  al.  Jl)  on  Rayleigh  wave. 

corrugation).  Rayleigh  and  Scholte  waves  propagating  along  an  isotropic  elastic  hatr-space  are  nondispersive.  The 
source/receiver  range  was  fixed  at  x  -  12  cm.  The  vertical  amplitude  scale  was  arbitrary.  Figure  1(b)  compares 
the  Rayleigh  wave  detected  at  two  different  locations  separated  by  one  cycle  along  the  corrugation  (x  ■=*  JO  cm  and 
x  **  12  cm).  The  lime  delay  between  the  two  waveforms  caused  the  corresponding  Bragg  frequency  component 
to  be  shifted  by  180°  (arrow  pointing  to  out  of  phase  components).  The  measured  Rayleigh  wave  Bragg  frequency 
along  this  sinusoidal  surface  was  2 1 .62  KHz.  The  periodic  surface  decreased  the  Rayleigh  wave  velocity  from  1 220 
m/s  to  865  m/s  at  the  Bragg  frequency.  The  corresponding  Scholte  wave  experimental  results,  measured  al  the  same 
two  locations,  are  shown  in  Figure  1(c).  The  Scholte  wave  . Bragg  frequency  was  20.2  KHz  indicating  that  the 
Scholte  wave  velocity  decreased  from  1 109  m/s  to  808  m/s. 

The  Scholte  and  Rayleigh  wave  experimental  dispersion  results  obtained  from  a  "soft '  elastic  models  are  compared 
with  replotted  scaled  numerical  results  1 1 )  on  Rayleigh  wave  dispersion  for  a  ’hard’  corrugated  solid  half-space. 
The  scaling  was  achieved  hy  simply  multiplying  the  dispersion  curve  by  the  Rayleigh  wave  velocity  ratio  between 
the  theoretical  "hard*  model  and  the  experimental  "soft*  model.  This  yields  approximate  results  because  Poisson's 
ratio  was  not  changed.  Figure  2  compares  the  Rayleigh  and  Scholte  wave  experimental  dispersion  results  with 
scaled  theoretical  dispersion  curve  of  Glass  et  al.  The  experimental  results  generally  agreed  with  the  Rayleigh  wave 
theoretical  results  by  Glass  et  al..  Work  is  in  progress  to  compare  the  experimental  results  with  Bostrom's 
calculated  results  [2]  on  Rayleigh  waves.  To  the  author's  knowledge,  there  are  no  known  published  theoretical  or 
numerical  results  characterizing  the  dispersion  of  Scholte  waves  by  a  large  amplitude  corrugated  liquid/solid 
interlace.  Quantitative  controlled  laboratory  ultrasonic  measurements  revealed  that  a  rippled  "soft*  liquid/solid 
interface  can  decrease  the  velocity  of  high-frequency  Scholte  wave  components  propagating  normal  to  the  ripples 
by  more  than  70%.  The  findings  indicate  the  potential  importance  of  seafloor  topography  on  acoustic  inversion  and 
detection  of  buried  objects  in  shallow  water. 
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4.2  Rayleigh  Waveform  from  Rippled  Drained  Wet  Sand  Surface 

The  effect  of  a  rippled  wet  sand  surface  on  Rayleigh  wave  dispersion  and  attenuation  is 
shown  in  Fig.  33.  The  top  trace  was  obtained  from  a  flat  surface  of  a  drained  sand  model. 
The  water  level  was  9  cm  below  the  sand  surface.  A  compressional  source  was  placed  on 
the  free  wet  sand  surface  to  generate  a  Rayleigh  wave.  The  receiver  was  a  2.5  cm 
diameter  shear  transducer  horizontally  oriented  and  buried  flush  with  the  sand  surface  to 
detect  both  shear  waves  and  Rayleigh  waves.  The  source-receiver  range  was  10  cm.  The 
source  was  excited  with  a  broadband  pulse  and  generated  a  transient  Rayleigh  wave  with 
a  3.1  KHz  center  frequency.  The  shear  wave  arrived  before  the  Rayleigh  wave  and 
partially  overlapped  the  Rayleigh  wave.  The  shear  wave  center  frequency  was  5  KHz,  and 
the  Rayleigh  wave  velocity  was  65  m/s.  A  sinusoidal  corrugation  was  made  on  the  wet 
sand  surface  with  a  1  cm  spatial  wavelength  and  0.25  cm  amplitude.  The  sand  ripples 
attenuated  the  Rayleigh  wave  and  generated  a  dispersive  wavetrain  with  a  2.32  KHz 
Bragg  frequency  [19].  The  sand  ripples  decreased  the  effective  Rayleigh  wave  velocity. 


Source 


Receiver 


Wet  Sand 


Rayleigh  Wave  on  Flat  Wet  Sand 


Rayleigh  Wave  on  Rippled  Wet  Sand 


Fig.  33.  Effect  of  sand  ripples  on  Rayleigh  wave  dispersion.  Corrugation  spatial 
wavelength  1cm,  ripple  amplitude  0.25  cm,  source-receiver  distance  10  cm.  Bragg 
frequency  resonance  2.32  KHz. 
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SECTION  5 


5.  Wedge  Waves  Along  Sand  Dollar  Skeleton  Edge 

In  order  to  reliably  detect  and  identify  small  buried  objects  in  littoral  sandy  sediments 
using  acoustic  waves,  it  is  essential  to  understand  the  seismoacoustic  response  of  partially 
buried  inclusions  such  as  seashells  and  sand  dollars.  This  section  presents  the  first 
experimental  results  on  wedge  waves  propagating  along  the  edge  of  sand  dollar  skeletons. 
Ash  et  al.  [34]  suggested  the  presence  of  a  slow  edge  wave  along  the  tip  of  a  triangular 
ridge.  In  1972,  Lagasse  [35]  discovered  a  nondispersive  flexural  edge  wave  propagating 
along  the  apex  of  a  small  angle  elastic  wedge  using  a  numerical  finite  element  approach. 
The  same  year,  Maradudin  et  al.  [36  ]  independently  presented  theoretical  results  on  edge 
modes  in  finite  crystals.  The  wedge  wave  is  a  unique  nondispersive  flexural  wave  that 
propagates  with  a  velocity  that  depends  on  the  wedge  apex  angle  according  to  Lagasse's 
empirical  relation 

Vm  =  Cr  sin  m0 

where  Vm  =  antisymmetric  mode  velocity,  Cr  =  Rayleigh  wave  velocity,  m  =  mode 
number  and  0  =  wedge  angle.  The  number  of  antisymmetric  edge  modes  is  determined  by 
the  relation  m0  <  90°.  These  waves  will  be  referred  to  here  as  "wedge  waves"  in  order 
not  to  confuse  them  with  the  dispersive  flexural  edge  waves  along  the  edge  of  a  plate. 
Chamuel  [37]  reported  on  edge  waves  along  immersed  elastic  elliptical  wedge  with  range- 
dependent  apex  angle.  The  sand  dollar  has  the  geometry  of  a  wedge  with  a  variable  angle. 
A  sand  dollar  skeleton  is  formed  from  two  fused  plates.  Fig.  34  shows  a  sand  dollar  cross 
section  revealing  the  wedge  geometry  near  its  edge.  The  wedge  angle  between  the  fused 
plates  varies  between  1 5 0  - 1 8 0 .  The  edge  is  also  rounded  with  a  1  -4  mm  truncation. 


Fig.  34.  Sand  dollar  skeleton  cross  section  showing  wedge  geometry  with  variable  edge  angle. 
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The  thickness  of  the  sand  dollar  edge  varies  along  the  perimeter.  Waveforms  of  wedge 
waves  obtained  frorti  the  thin  and  thick  edge  regions  are  compared  in  Fig.  35.  Two 
superimposed  sets  of  wedge  waves  are  shown  in  the  figure  demonstrating  the  close 
similar  elastic  properties  of  sand  dollars  of  the  same  age.  This  may  be  acoustically 
significant  in  single-age  group  sand  dollar  beds.  The  wedge  wave  velocity  is  greater  along 
the  thick  edge  portion 

Waveforms  from  a  5.04  cm  edge  segment  cut  from  a  sand  dollar  skeleton  with  a  diameter  of 
11.4  cm  are  compared  in  Fig.  36.  The  top  trace  was  obtained  with  the  source  and  receiver 
transducers  oriented  nearly  parallel  to  the  sand  dollar  edge.  A  longitudinal  wave  and  a  small 
edge  wave  were  observed.  The  middle  trace  shows  a  symmetric  Rayleigh  edge  wave  obtained 
with  the  transducers  oriented  normal  to  the  edge  and  parallel  to  the  wedge  center  plane.  The 
bottom  trace  shows  a  large  antisymmetric  flexural  wedge  wave  generated  and  detected  with 
the  transducers  normal  to  both  the  edge  and  center  plane.  The  measured  compressional, 
Rayleigh,  and  antisymmetric  wedge  waves  velocities  were  3702  m/s,  23 1 1  m/s,  and  1774  m/s, 
respectively. 

Waveforms  from  whole  uncut  sand  dollar  skeletons  are  presented  in  Fig.  37.  The  wedge  wave 
was  detected  at  one  location  along  the  edge  (top  trace),  and  then  the  receiver  was  moved  9  cm 
further  along  the  edge  (bottom  trace).  The  measured  Rayleigh  wave  velocity  was  2160  m/s 
and  the  wedge  wave  velocity  varied  along  the  wavepath  from  1572  m/s  to  1647  m/s.  Water 
loading  decreased  the  wedge  wave  velocity  and  amplitude.  Fig.  38  compares  wedge  waves 
from  a  sand  dollar  skeleton  in  air  and  in  water.  Further  research  is  needed  to  characterize  the 
seismoacoustic  response  of  seashells  and  sand  dollars.'. 


Thick  Section 


Thin  Section 


Fig.  35.  Superimposed  antisymmetric  edge  waves  from  two  sand  dollars  with  same  size 
showing  close  similarity  between  their  elastic  properties. 
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Fig.  36,  Comparison  of  waveforms  from  a  sand  dollar  edge  segment  5,04  cm  in  length. 
Top:  Longitudinal  edge  excitation  .  Middle:  Normal  symmetric  edge  excitation. 
Bottom:  Antisymmetric  flexural  edge  excitation. 


Fig.  37.  Wedge  wave  detected  at  two  locations  separated  by  9  cm  along  edge. 


Fig.  38.  Comparison  of  sand  dollar  skeleton  wedge  waves  in  air  and  in  water. 
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CONCLUSIONS 

A  summary  of  the  new  findings  was  given  at  the  beginning  of  the  report.  The  results 
revealed  gaps  in  fundamental  physical  understanding  of  littoral  seismoacoustic  phenomena 
in  sandy  sediments.  It  became  evident  that  existing  theories  do  not  predict  the  measured 
waveforms  in  water-covered  or  drained  sand.  Further  quantitative  research  is  needed  to 
characterize  the  complex  seismoacoustic  properties  of  various  sand  conditions  (water- 
saturated,  liquefied,  drained,  disturbed,  compacted,  etc...).  Some  of  the  questions  that 
need  to  be  answered  include;  What  are  the  characteristics  of  different  types  of 
seismoacoustic  waves  existing  in  water-covered  sand  and  partially-saturated  sand?  What 
are  the  effects  of  water  content  and  air  bubbles  in  sandy  sediments  on  broadband  pulsed 
compressiona!  waves?  What  are  the  seismoacoustic  properties  of  natural  sand  ripples 
compared  to  plowed  ripples?  What  is  the  compressional  wave  velocity  between  10  KHz 
and  70  KHz  in  sand  near  the  beach?  How  does  the  frequency  dependence  of  the 
compressional  wave  critical  angle  affect  ambient  noise  and  surfseism?  How  to  determine 
grain  to  grain  contact  conditions  and  measure  air  bubbles  content  from  broadband 
compressional  wave  data?  How  to  use  broadband  compressiona)  waves  to  monitor 
changes  in  sand  grain  to  grain  contact  while  coring?  What  are  the  effects  of  puddles  near 
the  beach  on  converting  Rayleigh  and  Scholte  waves?  What  are  the  effects  of  localized 
liquefied  regions  in  shallow  water  along  the  seafloor  on  transmission  loss?  What  is  the 
combined  dispersion  effect  of  ripples,  water  level,  and  depth-dependence  of  velocities? 
Should  the  resonance  frequency  of  bubbles  in  sediments  decrease  or  increase  compared  to 
bubbles  in  water?  Velea  et  al.  [17]  reported  that  for  a  2.5  mm  diameter  bubble,  the 
resonance  frequency  decreases  from  1300  to  800  Hz  in  a  natural  sediment  with  a  mean 
particle  size  of  30  nvn,  while  it  is  mentioned  in  Reference  [33]  that  in  general,  "the 
resonance  frequency  of  the  bubble  in  sediment  is  greater  than  for  the  equivalent  bubble  in 
water"  due  to  the  sediment  shear  modulus  term.  It  is  hoped  that  the  outcome  of  this 
research  would  lead  to  future  quantitative  studies  to  gain  fundamental  understanding  of 
acoustic  oceanography  and  achieve  reliable  seismoacoustic  detection  of  buried  objects 
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CONCLUSIONS 


A  summary  of  the  new  findings  was  given  at  the  beginning  of  the  report.  The  results 
revealed  gaps  in  fundamental  physical  understanding  of  littoral  seismoacoustic  phenomena 
in  sandy  sediments.  It  became  evident  that  existing  theories  do  not  predict  the  measured 
waveforms  in  water-covered  or  drained  sand.  Further  quantitative  research  is  needed  to 
characterize  the  complex  seismoacoustic  properties  of  various  sand  conditions  (water- 
saturated,  liquefied,  drained,  disturbed,  compacted,  etc...).  Some  of  the  questions  that 
need  to  be  answered  include:  What  are  the  characteristics  of  different  types  of 
seismoacoustic  waves  existing  in  water-covered  sand  and  partially-saturated  sand?  What 
are  the  effects  of  water  content  and  air  bubbles  in  sandy  sediments  on  broadband  pulsed 
compressional  waves?  What  are  the  seismoacoustic  properties  of  natural  sand  ripples 
compared  to  plowed  ripples?  What  is  the  compressional  wave  velocity  between  10  KHz 
and  70  KHz  in  sand  near  the  beach?  How  does  the  frequency  dependence  of  the 
compressional  wave  critical  angle  affect  ambient  noise  and  surfseism?  How  to  determine 
grain  to  grain  contact  conditions  and  measure  air  bubbles  content  from  broadband 
compressional  wave  data?  How  to  use  broadband  compressional  waves  to  monitor 
changes  in  sand  grain  to  grain  contact  while  coring?  What  are  the  effects  of  puddles  near 
the  beach'  on  converting  Rayleigh  and  Scholte  waves?  What  are  the  effects  of  localized 
liquefied  regions  in  shallow  water  along  the  seafloor  on  transmission  loss?  What  is  the 
combined  dispersion  effect  of  ripples,  water  level,  and  depth-dependence  of  velocities? 
Should  the  resonance  frequency  of  bubbles  in  sediments  decrease  or  increase  compared  to 
bubbles  in  water?  Velea  et  al.  [17]  reported  that  for  a  2.5  mm  diameter  bubble,  the 
resonance  frequency  decreases  from  1300  to  800  Hz  in  a  natural  sediment  with  a  mean 
particle  size  of  30  pm,  while  it  is  mentioned  in  Reference  [33]  that  in  general,  "the 
resonance  frequency  of  the  bubble  in  sediment  is  greater  than  for  the  equivalent  bubble  in 
water"  due  to  the  sediment  shear  modulus  term.  It  is  hoped  that  the  outcome  of  this 
research  would  lead  to  future  quantitative  studies  to  gain  fundamental  understanding  of 
acoustic  oceanography  and  achieve  reliable  seismoacoustic  detection  of  buried  objects. 
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Scholte  Wave  Dispersion  by  Rippled  Liquid/Solid  Interface  Topography 

Jacques  R.  Chamuel 
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Abstract:  Quantitative  experimental  results  are  presented  characterizing  the  dispersion  of  Scholte  waves  propagating  along 
a  rippled  surface  of  an  immersed  "soft"  solid  half-space.  One  model  has  a  sinusoidal  profile  with  a  ratio  of  corrugation 
amplitude  to  wavelength  of  0.25.  Scholte  and  Rayleigh  wave  experimental  dispersion  results  are  compared  with  numerical 
results  on  Rayleigh  wave  dispersion  along  a  “hard"  corrugated  solid  half-space  [N.E.  Glass,  R,  Loudon,  and  A.  A. 
Maradudin,  Physical  Review  B  24,  6843-6861  (1981)].  The  findings  reveal  that  a  rippled  "soft"  liquid/solid  interface  can 
decrease  the  velocity  of  high-frequency  Scholte  wave  components  propagating  normal  to  the  ripples  by  more  than  70%. 


INTRODUCTION 

The  effect  of  a  rippled  liquid/solid  interface  on  the  dispersion  of  Scholte  waves  is  not  known.  Frequency- 
dependent  multiple  scattering  phenomena  occur  as  near-grazing  underwater  acoustic  waves  and  Scholte  waves 
interact  with  periodic  seafloor  topography  having  dimensions  comparable  to  the  wavelength.  In  order  to  develop 
physical  understanding  of  the  interaction  of  broadband  transient  underwater  acoustic  waves  and  Scholte  waves  with 
a  rippled  ocean  bottom,  experimental  ultrasonic  experiments  were  conducted  on  immersed  isotropic  "soft"  models 
with  a  sinusoidally  corrugated  surface.  Quantitative  experimental  results  were  obtained  characterizing  the  dispersion 
of  Rayleigh  and  Scholte  waves  caused  by  a  periodic  iiquid/solid  interface. 


MODEL  EXPERIMENTS  AND  RESULTS 


A  sinuoidal  pattern  was  machined  along  the  edge  of  a  4.5  mm  thick  Plexiglas  plate.  The  spatial  wavelength  of 
the  periodic  pattern  was  2  cm  and  the  peak-to-peak  amplitude  was  1  cm.  The  measured  Plexiglas  velocities  were 
as  follows:  compressional  plate  velocity  =  2313  m/s,  shear  wave  velocity  =  1344  m/s,  flat  edge  Rayleigh  wave 
velocity  =  1220  m/s,  flat  edge  Scholte  wave  velocity  =  1 109  m/s.  The  top  trace  of  Figure  1(a)  compares  the 
dispersed  waveforms  of  Rayleigh  and  Scholte  waves  propagating  along  the  periodic  surface.  Rayleigh  and  Scholte 
waves  were  obtained  from  air/Plexiglas  and  water/Plexiglas  models,  respectively.  The  Scholte  waveform  was 
stretched  more  to  the  right  than  the  Rayleigh  wave  indicating  that  it  was  slower  than  the  Rayleigh  wave.  The 
bottom  trace  of  Figure  1(a)  shows  the  arrival  time  of  the  corresponding  Rayleigh  wave  along  a  flat  surface  (no 


(a)  (b)  (c) 


FIGURE  1.  Comparison  of  received  waveforms  and  arrival  time;  (a)  Top  Trace:  Rayleigh  and  Scholte  waves 
along  corrugated  surface  (range  x  =  12  cm),  Bottom  Trace:  reference  Rayleigh  wave  from  flat  solid  surface,  (b) 
Rayleigh  wave  along  cormgated  surface  detected  at  two  positions  separated  by  one  corrugation  cycle  (2  cm),  (c) 
Scholte  wave  along  corrugated  surface  detected  at  the  same  two  positions  in  (b). 
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Experimental  [Chamuel  1997) 


Rayleigh  Wave 
Rayleigh  Wave 
Scholte  Wave 


FIGURE  2.  Rayleigh  and  Scholte  wave  dispersion  experimental  results  compared  with  approximately  scaled 
theoretical  dispersion  results  by  Glass  et  al,  [1]  on  Rayleigh  wave. 


corrugation).  Rayleigh  and  Scholte  waves  propagating  along  an  isotropic  elastic  half-space  are  nondispersive.  The 
source/receiver  range  was  fixed  at  x  =  12  cm.  The  vertical  amplitude  scale  was  arbitrary.  Figure  1(b)  compares 
the  Rayleigh  wave  detected  at  two  different  locations  separated  by  one  cycle  along  the  corrugation  (x  =  10  cm  and 
X  =  12  cm).  The  time  delay  between  the  two  waveforms  caused  the  corresponding  Bragg  frequency  component 
to  be  shifted  by  180°  (arrow  pointing  to  out  of  phase  components).  The  measured  Rayleigh  wave  Bragg  frequency 
along  this  sinusoidal  surface  was  21.62  KHz.  The  periodic  surface  decreased  the  Rayleigh  wave  velocity  from  1220 
m/s  to  865  m/s  at  the  Bragg  frequency.  The  corresponding  Scholte  wave  experimental  results,  measured  at  the  same 
two  locations,  are  shown  in  Figure  1(c).  The  Scholte  wave  Bragg  frequency  was  20.2  KHz  indicating  that  the 
Scholte  wave  velocity  decreased  from  1109  m/s  to  808  m/s. 

The  Scholte  and  Rayleigh  wave  experimental  dispersion  results  obtained  from  a  "soft"  elastic  models  are  compared 
with  replotted  scaled  numerical  results  [1]  on  Rayleigh  wave  dispersion  for  a  "hard"  corrugated  solid  half-space. 
The  scaling  was  achieved  by  simply  multiplying  the  dispersion  curve  by  the  Rayleigh  wave  velocity  ratio  between 
the  theoretical  "hard"  model  and  the  experimental  "soft"  model.  This  yields  approximate  results  because  Poisson’s 
ratio  was  not  changed.  Figure  2  compares  the  Rayleigh  and  Scholte  wave  experimental  dispersion  results  with 
scaled  theoretical  dispersion  curve  of  Glass  et  al.  The  experimental  results  generally  agreed  with  the  Rayleigh  wave 
theoretical  results  by  Glass  et  al..  Work  is  in  progress  to  compare  the  experimental  results  with  Bostrom’s 
calculated  results  [2]  on  Rayleigh  waves.  To  the  author’s  knowledge,  there  are  no  known  published  theoretical  or 
numerical  results  characterizing  the  dispersion  of  Scholte  waves  by  a  large  amplitude  corrugated  liquid/solid 
interface.  Quantitative  controlled  laboratory  ultrasonic  measurements  revealed  that  a  rippled  "soft"  liquid/solid 
interface  can  decrease  the  velocity  of  high-frequency  Scholte  wave  components  propagating  normal  to  the  ripples 
by  more  than  70%.  The  findings  indicate  the  potential  importance  of  seafloor  topography  on  acoustic  inversion  and 
detection  of  buried  objects  in  shallow  water. 
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